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1 

COMPUTER-VISION SYSTEM FOR CLASSIFICATION AND SPATIAL LOCALIZATION OF 
BOUNDED 3D-OBJECTS. 

5 

FIELD OF THE INVENTION 

The invention relates to a method for object recognition in a computer vision system, more 
specifically the method relates to classification and spatial localization of bounded 3D- 
objects. 

10 

BACKGROUND OF THE INVENTION 

A bottleneck in the automation of production processes is the feeding of components and 
semi-manufactured articles to automatic systems for machining, assembly, painting, 
packing etc. Three main types of systems are available today: 1) vibration bowls, 2) 

15 fixtures, and 3) computer vision systems. Vibrating bowls are suitable only for components 
of small dimensions (less than about 5 cm). Fixtures are expensive, since the entire 
internal storage must be based on such fixtures. Both types of systems must be 
redesigned and remanufactured when new components are introduced. The computer 
vision systems developed so far have serious drawbacks. Some systems have 

20 unacceptably low processing speeds, others have poor generality. The fast and general 
systems available today require the objects to lie scattered on a flat conveyer belt, and the 
object-camera distance must be much larger than the object height. The latter limitation is 
fundamental for the present systems, as the recognition model used does not include 
perspective effects in the 3D-2D transformation of the camera. Thus, for parts higher than 

25 5 - 10 cm, standard computer vision systems demand inconveniently remote cameras. 
Furthermore, they are not able to guide robots to structured grasping randomly oriented 
parts piled in boxes and pallets. 

Another bottleneck is present when recycled articles are to be classified as they arrive to 
30 the recycling plants. The rebuilding of parts used in consumer products, particularly in 

cars, is expected to increase in the future for environmental and resource reasons. Prior to 
the rebuilding process there is a need for classification. 

A third example of a field with insufficient technology at present is fast navigation of 
35 mobile robots in structured environments. The camera based navigation systems require 
recognition of building elements, stationary furniture etc. Segments of these can be 
considered to be bounded 3D objects. 
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Furthermore the system can be used in satellite applications for identification and 
classification of vehicles, buildings etc. 

SUMMARY OF THE INVENTION 

5 

According to preferred embodiment of the invention, recognition and/or localisation of 
objects is based on primitives identified in a recognition image of an object. Thus, in a first 
aspect, the present invention relates to a method of determining contours, preferably level 
contours and primitives in a digital image, said method comprising the steps of: 
10 - generating the gradients of the digital image; 

finding one or more local maxima of the absolute gradients; 

use the one or more local maxima as seeds for generating contours, the generation 
of the contours for each seed comprising determining an ordered list of points 
representing positions in the digital image and belonging to a contour; 
15 - for all of said positions determining the curvature, preferably determined as dO/ds 
preferably pixel units, of the contours; 

from the determined curvatures determine primitives as characteristic points on or 
segments of the contours. 

20 Based on the primitives derived from training image recognition and/or localisation of an 
object may preferably be performed by a method according to a second aspect of the 
present invention, which second aspect relates a method of recognition, such as 
classification and/or localisation of three dimensional objects, said one or more objects 
being imaged so as to provide a recognition image being a two dimensional digital image 
25 of the object, said method utilises a database in which numerical descriptors are stored for 
a number of training images, the numerical descriptors are the intrinsic and extrinsic 
properties of a feature, said method comprising: 

identifying features, being predefined sets of primitives, for the Image 
extracting numerical descriptors of the features, said numerical descriptors being of 
30 the two kind: 

extrinsic properties of the feature, that is the location and orientation of the 
feature in the image, and 

intrinsic properties of the feature being derived after a homographic 
transformation being applied to the feature 
35 - matching said properties with those stored in the database and in case a match is 

found assign the object corresponding to the properties matched in the database to 
be similar to the object of the object to be recognised. 
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In a third aspect the present invention relates to a method of generating a database useful 
in connection with localising and/or classifying a three dimensional object, said object 
being imaged so as to provide a two dimensional digital image of the object, 
said method utilises the method according to the first and/or the second aspect of the 
5 invention for determining primitives in the two dimensional digital image of the object, said 
method comprising: 

identifying features, being predefined sets of primitives, in a number of digital 
images of one or more object, the images represent different localisations of the 
one or more object; 

10 - extracting and storing in the database, numerical descriptors of the features, said 

numerical descriptors being of the two kind: 

extrinsic properties of the feature, that is the location and orientation of the 
feature in the image, and 

intrinsic properties of the feature being derived after a homographic 
15 transformation being applied to the feature. 

The present invention thus allows for the system to recognize, classify and localize objects. 

The invention may furthermore comprise the step of eliminating potential seed points 
20 identified near already defined contours. This is done preferably in order to avoid 
generation of contours that are too close to already existing contours. 

Furthermore the generation of the contours may comprise assigning the list of points 
representing positions in the digital image, each point having a value being assigned to be 
25 common with the value of the seed. 

Even further the generation of contours may be defined as determining an ordered list of 
points comprising points following in each point the direction of the maximum gradient. 

30 Even further the generation of the contours may comprise assigning the list of points 
following in each point the direction of the maximum or minimal gradient detected 
perpendicular to a contour direction. Which gradient to follow may be decided upon which 
contour direction that is chosen. 



35 



Moreover the generation of the contours may comprise assigning a list of pixels with 
values being above or below the value of the seed and one or more neighbour pixels with 
value below or above said value of the seed. 
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The list of pixels Is preferably established by moving through the digital image in a 
predetermined manner. The established list may be an ordered list of pixels which would 
enhance the speed of searching the list since a certain value such as a max or min would 
e.g. be in the top respectively in the bottom of the list. However other solutions may also 
5 be applied. 

Moreover the contours may be determined from an interpolation based on the list of pixels. 
This is preferably done in order to obtain a smoother contour. 

10 The creation of a gradient picture may be achieved by determining gradients by calculating 
the difference between numerical values assigned to neighbouring pixels. In this way a 
gradient picture is obtained from which further information may be extracted. 

The gradients may be stored in an array in which each element preferably corresponds to a 
15 specific position in the first image and being a numerical value representing the value of 
the gradient of the first image's tones in the specific position. 

The curvatures which may be used for generating primitives are preferably established as 
K=de/ds where 9 is the tangent direction at a point on a contour and s is the arc length 
20 measured from a reference point. 

The primitives mentioned in the first, second and third aspect above preferably comprise of 
one or more of the following characteristics: 

segments of straight lines, 
25 - segments of relatively large radius circles, 

inflection points, 

points of maximum numerical value of the curvature, said points being preferably 
assigned to be corners, 

points separating portions of very low and very high numerical value of the 
30 curvature, and 

small area entities enclosed by a contour. 

The generated contours mentioned above may be searched for one or more of the 
following primitives: 

35 - inflection point, being a region of or a point on the contour having values of the 
absolute value of the curvature being higher than a predefined level; 
concave corner, being a region of or a point on the contour having positive peaks of 
curvature; 
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convex corner, being a region of or a point on the contour having negative peaks of 
curvature; 

straight segment, being segments of the contour having zero curvature; 

and/or 

5 - circular segments, being segments of the contour having constant curvature. 

According to the second aspect of the invention which relates to a method of recognition, 
the matching step may further comprise matching a recognition image with training 
images stored in a database, the method for matching may comprise the following steps: 
10 - for each training image: 

determining the values of roll, tilt and pan of the transformations bringing 
the features of the recognition image to be identical with the features of the 
training image; 

identify clusters in the parameter space defined by the values of roll, tilt 
15 and pan determined by said transformations 

and 

identify clusters having predefined intensity as corresponding to an object type and 
localisation. 

20 The database may comprise for each image one or more records each representing a 
feature with its intrinsic properties and its extrinsic properties. 

Furthermore the matching may comprise the steps of: 
resetting the roll, tilt and pan parameter space, 
25 - for each feature in the recognition image, matching properties of the recognition 
image with the properties stored in the database, 

in case of match: determining roll, tilt, and pan based on the extrinsic 
properties from the database and from the recognition image, 
updating the parameter space, and 
30 - test for clustering and store coordinates of clusters with sufficiently high 

density/population with an index of the training image, 
repeating the steps until all features in the recognition image have been matched. 

The determination of the roll, tilt and pan are preferably only done for features having 
35 similar or identical intrinsic properties compared to the intrinsic properties in the database. 
In this way a lot of computing time may be saved. 
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Moreover the matching may comprise comparing the intrinsic descriptors of the recognition 
Image with the intrinsic descriptors stored in the database thereby selecting matching 
features. 

5 The generation of database mentioned in the third aspect may furthermore be generated 
according to the steps and features described below. 

Preferably the extrinsic properties comprises a reference point and a reference direction. 
Whereas the intrinsic properties preferably comprises numerical quantities of features. 

10 

The object may be imaged by a method using at least two imaging devices thereby 
generating at least two recognition images of the object and wherein the method for 
matching may be applied to each recognition image and wherein the match found for each 
recognition image are compared. 

15 

Furthermore the method using two imaging devices may comprise the steps of: 

for each imaging device, providing an estimate for the three dimensional reference 
point of the object, 

for each imaging device, calculating a line from the imaging device pinhole to the 
20 estimated reference point, 

and when at least two or more lines have been provided, 

discarding the estimates in the case that the said two or more lines do not 
essentially intersect in three dimensions, 
and when the said two or more lines essentially intersect, 
25 - estimating a global position of the reference point based on the pseudo intersection 

between the lines obtained from each imaging device. 

In the following, the invention and in particular preferred embodiments thereof, will be 
presented in greater details in connection with the accompanying drawings. 

30 

BRIEF DESCRIPTION OF FIGURES 

Figure 1, illustrates a tilt-pan homographic transformation. 
35 Figure 2a - 2c, shows primitives, pairs of primitives and angles. 
Figure 3, shows an example of an image in four different windows. 
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Figure 4, shows the contours in the upper right corner of the image in Figure 3. The 
window above in figure 4 shows the contour with subpixel accuracy while the window 
below in figure 4 shows the integer pixel positions of the contour. 

5 Figure 5, illustrates Curvature k(s) (in radians/pixel) as a function of arc length s (in 

pixels) along one of the contours of Figure 3, Window 4. The symbols are used for showing 
the correspondence, see figure 6. 

Figure 6, illustrates the outer contour found for the image shown in Figure 3. The symbols 
10 used for characteristic features correspond to those given in Figure 5. 

Figure 7, illustrates a 3D triple brick object treated in the example for the pose 
determination. 

15 Figure 8, illustrates training photos of a Lego model of the triple brick structure. 

Figure 9, illustrates the curvature in units of radians/pixel as a function of the arc length 
(in pixels) along the contour in figure 8A. 

20 Figure 10, shows a combined view of a training image and a recognition image. 

Figure 11, shows a flow chart describing the processing for training. 

Figure 12, show a flow chart describing the processing for recognition. 

25 

Figure 13, illustrates two-camera operation. 
Figure 14, illustrates a pinhole model of a camera. 
30 Figure 15, illustrates the greytone landscape derived from an image. 
Figure 16, illustrates the training geometry. 

Figure 17, illustrates the structure of the database of descriptors derived from training 
35 images. 

Figure 18, illustrates the structure of descriptors derived from recognition image. 
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DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS OF THE INVENTION 

The invention described here is aimed at all the situations described in the section 
Background of the invention. Focus is put on the following properties: 
5 - Simple generation of training information 

Reasonably low volume of training information. 
Exact treatment of the perspective effects 
- Generality concerning the shape and the visual appearance of the objects, e.g. sharp 
3D edges and landmarks are not necessary 
10 - High speed recognition without extensive 2D matching between images or 3D 
reconstruction. 

Functionality 

The computer vision system is used for classifying and/or locating bounded 3D objects 
15 belonging to distinct classes. The system consists of one or more cameras whose images 
are interpreted in terms of 1) class of 3D objects and 2) their spatial position and 
orientation (pose). Its function is to some extent independent of possible partial occlusion 
by other objects and of poor image segmentation. The objects need not have characteristic 
decorations or sharp edges. The function is independent of camera position and object size 
20 relative to the object-camera distance. The image interpretation is speed optimized which 
implies the use of digital camera, swift electronic transmission of image data, and 
optimized code. Furthermore the camera used in the system does not necessarily have to 
be an optical camera, the camera can be of other kinds such as a thermal camera. 

25 Definitions 

Pinhole model of a camera, illustrated in figure 14: The frame (coordinate system) of the 
camera is defined by the axes u,v,w. The focal point has the coordinates (u,v,w)=(0,0,f), 
where f is the focal length of the camera. Preferably units of u,v and f are pixel units. 
Physical cameras have negative values of f. The relevant homographical transformation 
30 can be described by two successive rotations about the tilt-axis (parallel to the u-axis), 
and the pan axis (parallel to the v-axis). 

A camera is an imaging device with a pinhole, i.e. center of a perspective 3D-2D 
transformation, and an image plane. The optical axis is a line through the pinhole, 
35 essentially perpendicular to the image plane. The image of the optical axis of the camera is 
called the focal point, illustrated in figure 14. The image has two axes, a vertical (v) axis 
and a horizontal (u) axis. 
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Preferably the following 2D properties of the visual appearance of objects are considered: 

1) the outer contour (always existing), 

2) contours appearing inside the outer contour, 

3) images of sharp 3D edges of the object appearing inside the contour, and 
5 4) 2D edges in decorations. 

All these properties may appear as (one-dimensional) lines or curves in the image. In the 
following these features are called characteristic curves. Specific features of characteristic 
curves are called primitives. Primitives can be point-like (inflection points, points of 
10 maximum curvature etc.) or one-dimensional (straight sections, sections with constant 
curvature etc). Specific pairs, triplets, or higher sets of primitives are called features. The 
most useful types of features are pairs of primitives. A few of these are illustrated in figure 
2a and 2b. 

15 An image of a single, specific object and with a known object-camera pose taken by a 
specific camera is called a training view. An image of a scene to be interpreted by the 
system is called a recognition view. 

Numerical descriptors describes the intrinsic properties and extrinsic properties of a 
20 feature. Intrinsic properties are described by the rotation invariant descriptors of features, 
whereas the extrinsic properties are described by the location and rotation of a feature in 
an image. 

A feature preferably has three extrinsic descriptors: The two coordinates of the reference 
25 point of the feature, and the reference direction after homographic transformation. 

Level contours: A level contour is preferably an ordered list of image coordinates 
corresponding to a constant greytone value g. The coordinates are obtained by linear 
interpolation between two pixels, one with greytone above g, the other with greytone 
30 below g. 

Edge contours: An edge contour is preferably an ordered list of points representing 
positions in the digital image and belonging to a contour. A point is preferably obtained 
from a previous point by calculating the gradient perpendicular to the present motion 
35 direction and moving in the direction of the maximum gradient. Subpixel coordinates of 
contour points are obtaied by calculating the average greytone over (integer) pixel 
positions some distance behind and ahead of each contour point. This average greytone is 
used as a basis for an interpolation leading to a shifted position of the considered point. 
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Views/images or sections of views/images can be subject to 2D transformations. The 
transformations considered here are characterized by virtual rotations of the camera about 
its pin-hole. These transformations are denoted as homographic transformations. 
Homograpical transformations can be specified by successive camera rotations about 
5 specific axes. In a common notation tilt is a rotation about the axis parallel to the 

horizontal image axis, pan is a rotation about the axis parallel to the vertical image axis, 
and roll is a rotation about the optical axis. These rotations are illustrated in figure 14. Let 
n be an image or an image section. The transformed image or image section has the 
symbol & = H(n), where H is the homographic transformation. Any given point Q in the 

10 image defines a class of homographic transformations with the property that the point Q is 
transformed into the focal point. The image or image section after such a transformation 
has the symbol H Q (a). One member of this class of transformations H Q is characterized by 
a tilt followed by a pan, and no roll. This transformation will be called the tilt pan 
transformation H Qftp . There exist many other members of this class. It is preferred that 

15 they have well defined algorithms. 

In figure 1, a tilt-pan homographic transformation is illustrated. The original image before 
the transformation is the upper image. Below is the Image after the tilt-pan homographic 
transformation, wherein the tip of the dome is moved to the focal point. 

20 

Preferably the objects of the same class are uniform with respect to geometric form and to 
some extent also decoration and the non-occluded part of objects has sufficient 
characteristic curves. In order to achieve the best detection of characteristic curves, the 
illumination of the scene is preferably reasonably constant. 

25 

Overall description of the method 

The recognition is based on the analysis of a large number of training views. These training 
views are recorded by a camera viewing a real object or constructed using the CAD 
representation of the object. Characteristic curves are derived from the training views, and 
30 primitives of the curves are detected. Intrinsic and extrinsic descriptors of features are 
stored in a database together with data about the object class and pose of the view. The 
above activities related to training are performed off-line. 

A similar image analysis is performed during recognition. The remaining part of the 
35 recognition takes place in two stages: First the intrinsic descriptors of the recognition view 
are compared with those of the database. Second, among the best matching features it is 
explored which features agree mutually in the sense that they suggest the same object 
class at the same pose. 
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Methods for reduction of number of training views 

As a rigid body has 6 degrees of freedom, the diversity of views is very large. Two 
methods for reduction of the training volume are employed. First, the extrinsic descriptors 
5 are derived from tilt-pan homographically transformed images. The transformation used 
for a given feature is H Qftpf where Q is a reference point of the actual feature. Second, the 
intrinsic descriptors used in the match search are invariant to image rotations, equivalent 
to rolling the camera. The above two strategies imply that the volume of training views can 
be limited to three degrees of freedom in spite of the fact that a rigid object has six 
10 degrees of freedom. For each feature the training database preferably contains: 

a) Descriptors invariant to tilt-pan homographic transformation and to a roll operation, 

b) A rotation-describing descriptor for the angular 2D orientation relative to an image 
axis, 

15 c) The tilt and pan angle involved in the tilt-pan homographic transformation. 

Point a) requires that a reference point can be assigned to the feature. 
Point b) requires that a reference direction can be assigned to the feature. 

20 In the training session the reference direction and reference point may be assigned 
manually by a user. 

The three degrees of freedom involved in the training can be chosen to be the spherical 
pinhole coordinates (p, cp, 9) in the object frame, see figure 16. During training the optical 

25 axis is going through the origin of the object frame, and the roll angle of the camera is 
zero. Thus, p, cp, 9 are length, azimutal angle and horizontal angle, respectively, of the 
vector from the origin of the object frame and the pinhole. A user assigns the origin of the 
object frame in the training step. The intervals and step sizes of p, <p, and 9, necessary to 
be trained depend on the application. In case of moderate or weak perspective only few 

30 values of p needs to be trained as the linear dimensions of the features are approximately 
inversely proportional to p. 

Recognition: Transformation, match search, backtransformation and cluster analysis. 

In the first step of the recognition the recognition view is analyzed, the descriptors of 
35 transformed features are derived, and an appropriate number of the best matches 
between descriptors from the recognition view and those of the database (item a) in 
Section {Methods for reduction of number of training views) are found. In the second step 
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one considers the items b) and c), Section (Methods for reduction of number of training 
Wews), belonging to the recognition view and the matching record of the database. These 
data are used in a suitable backtransformation thereby calculating the full 3D pose 
suggested by the actual features. Clusters of suggestions (votes) in the 6-dimensional 
5 configuration space (one space for each object class) are interpreted as real objects. This 
cluster analysis is essential for eliminating false features, i.e. detected combination of 
primitives belonging to different objects. 

Primitives of characteristic curves and features 

10 Examples of preferred primitives for recognition: 

1) Segments of straight lines 

2) Segments of relatively large radius circles 

3) Inflection points 

4) Points of maximum curvature 

15 5) Points separating portions of very low curvature and very high curvature. 
6) In case of curves enclosing small areas: The 2D center-of-mass of this area. 

Figure 2 - 2a shows primitives mentioned in Section (Primitives of characteristic curves 
and features) above. Figure 2a shows examples of primitives, figure 2b shows pairs of 
20 primitives, their reference points (thin circles) and their reference directions (arrows). 
Figure 2c illustrates angles wherein, the angles r are rotation invariant, and the angles d 
are rotation-describing descriptors. 

The sets of primitives used in the system should preferably have the following properties: 
25 - a reference point 

- a reference direction (without 180 degree ambiguity) 

- one or more rotation invariant descriptors suitable for a match search. 

Any combination of two or more primitives fulfilling these conditions can be employed. 
30 Figure 2b shows example of suitable pairs of primitives including reference points and 

reference directions. In the case that segments of straight lines or circles are involved in a 
feature, then the recognition allowing partial occlusion of the features should preferably 
involve appropriate inequalities. 

35 Rotation invariant descriptors of the pairs of primitives are for example distances between 
point like primitives, angles between portions of straight lines, angles between tangents 
and lines connecting point-like primitives, etc. Figure 2c shows examples of rotation 
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Invariant angles, and the rotation describing angle (item b) in section (Methods for 
reduction of number of training views). 

Advantages using two or more cameras 

5 An uncertain component of the pose In single camera applications is p, I. e. the distance 
between the pinhole and the reference point of the object. Errors come from pixel 
discretization, camera noise and fluctuating object dimensions. The uncertainty can be 
reduced significantly by correlating findings from two or more cameras as follows. Each 
camera gives an estimate for the 3D reference point of the object. With uncertain p each 
10 camera defines a 3D line of high probability for the reference point position. The pseudo 
intersection between such lines is the most probable position of the reference point of the 
object. This is illustrated in figure 13. 

This method is related to stereovision. Conventional stereovision has a fundamental 
15 limitation since a too short base line (distance between pinholes) gives an inaccurate depth 
determination, while a too large base line (and large angles between optical axes) makes 
the identification of corresponding points/features difficult. In the presently introduced 
method using features matching those of a muiti view database, there is no need for 
finding corresponding points in the views. Therefore, the depth estimation achieved by a 
20 muiti camera version of the present invention is more accurate than with ordinary 
stereovision. 

Another advantage obtained by using more than one camera is the elimination of 
misclassifications and wrong pose estimations. This elimination is particularly important in 
25 case of objects with a symmetry plane viewed with weak perspective. 

Estimate of training volume and recognition times 

In a typical application, the step size in <p-0 space is 4 degrees. This implies about 3000 
views per value of p for an unrestricted angular range. Most applications need only 3-4 

30 different p-values, giving a total of about 10,000 images. A typical number of sets of 
primitives in each training view is 50, and the typical number of 4 byte floating-point 
entities in each database record is 8. Then the total volume of the database is of the order 
16 MByte for one object class. A speed optimized match search in this database is 
expected to last less than one second per object class on a 1GHz CPU. In applications 

35 where it is known a priori, that the object pose is confined to a smaller part of the (p,cp,9)- 
space, the above numbers can be reduced correspondingly. 
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The embodiment described herein comprises the preferred steps used in the computer 
vision system. The system is able to classify and locate 3D objects lying at random in front 
of one or two computer-vision cameras. As outlined in the summary of invention of the 
system, the recognition is preferably based on: 

5 

Determination of characteristic curves in the training and recognition images 

- Derivation of feature descriptors (primitives and pairs of primitives), 

- A recognition processes to be used in the 3D interpretation 

10 The characteristic curves used are edges In greytone images. In this description the 
'edges' are defined as level contours (curves of constant greytone) or edge contours 
(curves of high gradient). The method for deriving level contours is described and 
exemplified in Sect. (Derivation of level contours from greytone images). By using subpixel- 
defined contours, it is possible to derive reliable characteristic contour primitives (straight 

15 segments, inflection points, corners, etc.) as outlined in Sect. (Derivation of primitives and 
features from contours). The 3D interpretation using primitives derived from training and 
recognition images is described and exemplified in Sect. (Steps in the 3D interpretation). 

Derivation of contours from greytone images 

This section describes the Image analysis leading to level and edge contours. The definition 
20 of level contours: The greytone landscape of the frame in the upper right section of figure 
15a is shown in figure 15b. A level contour is preferably an ordered list of image 
coordinates corresponding to a constant greytone value g. The coordinates are obtained by 
linear interpolation between two pixels, one with greytone above g, the other with 
greytone below g. Below follows explanations of a few definitions. 

25 

1. A greytone image consists of a 2D array G[x,y] of greytones. Each array member is a 
pixel. 

2. Each pixel has integer coordinates in the image plane. 

3. In analogy with a landscape the greytones are considered as heights in the greytone 
30 landscape, this is illustrated in figure 15a and 15b. 

4. With suitable interpolation the greytone can be considered as a function of continuous 
pixel coordinates. 

5. A curve in the image plane going through points with common greytone g is called a 
level contour. Note that level contours preferably do not cross each other. 

35 6. The 'gradient' at a point (x,y) is defined as; max ( |G[x,y-l]-G[x,y+l]| , |G[x-l,y]- 
G[x+l,y]|). 

7. Pieces of level contours with high gradient are x edge like'. 
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It is the aim of this section to describe an efficient way of deriving meaningful level 
contours. The result of the image analysis is a list of level contour segments, and each 
level contour segment is a list of pixel positions. 

5 Deriving seeds for contours 

In the first step a "gradient image' is derived as an image in which the greytone is equal to 
the gradient of the original image. A potential seed is defined as local maximum in the 
gradient image. A list of potential seeds is formed. Maxima with gradients below a 
threshold are not used as seeds. The list of seeds contains the greytone, the gradient and 
10 the pixel coordinates. This list is sorted according to gradient magnitude. Figure 3 shows 
an example of Image analysis leading to seeds. Window 1 in figure 3 shows the original 
image, Window 2 in figure 3 shows the gradient image. Window 3 in figure 3 shows the 
potential seeds. Window 4 in figure 3 shows the contours derived. 

15 Deriving the contours 

The first level contour to be generated uses the seed with highest gradient. A standard 
contour search is applied by using the greytone threshold equal to the greytone of the 
seed. The contour is followed until at least one of the following stop conditions is fulfilled: 
1) the image border is reached, 
20 2) the seed is reached again (closed contour) or 

3) the absolute value of the gradient at the next contour point falls below a threshold. 

The contour search is preferably bi-directional unless the contour is closed. Potential seeds 
closer than 1-2 pixels to the contour derived are disabled/eliminated. 

25 

The pixel positions in all contours may be shifted according to a linear interpolation using 
the greytone value characteristic for each contour. The result is shown in figure 4. 

The windows in figure 4 shows the contours in the upper right comer of the image of 
30 Figure 3. The window above in figure 4 shows the contour with subpixel accuracy while the 
window below in figure 4 shows the integer pixel positions of the contour. 

The next contour is generated using the non-disabled seed with highest gradient. New 
contours are then repeatedly generated until the list of seeds is exhausted. Figure 3, 
35 Window 4, shows an example of contours derived. In this example, the number of contour 
sections is 9. Weaker contours than those shown, can be generated by choosing a smaller 
value for the minimum gradient of seeds. 
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Preferably the following procedure and constraints are applied: 

a) The level contours are drawn starting from a seed point. 

b) Potential seeds are pixels with local maximum of the gradient magnitude. 

c) The level contours are derived in a succession of decreasing gradient magnitude of 
5 their seeds. Starting with the seed having the highest gradient magnitude. 

d) Level contours positions are initially found as integer pixel positions with greytones 
above the value g and at least one neighbour pixel with greytone below the value g. 
Interpolated pixel positions are obtained by shifting each (integer) pixel position to new 
position derived by interpolation. The position list in each level contour is ordered so 

10 that neighbour indices in the list correspond to neighbour positions in the image. 

e) When moving along a direction with increasing position index, then regions with 
greytones higher than g are preferably at the right hand side. 

In a second embodiment so-called edge contours are generated. In this case the above 
15 procedure and constraints may be applied with a change in step d, the new step d may be 
performed as described below: 

Edge contours positions are found as pixel positions starting from a seed (first pixel) and 
moving towards the next pixel having the maximum gradient perpendicular to the present 
20 contour direction, and repeating these steps until one of the stop conditions is met. Sub 
pixel positions of edge contour points are preferably derived by calculating the average 
greytone over some pixels behind and ahead of each contour position. This average 
greytone may be used as a basis for interpolation. 

Derivation of primitives and features from contours 

25 A primitive is a point on, or a segment of a contour with characteristic behaviour of the 
curvature, see figure 2a. The primitives listed in the summary of the invention are: 

a) Segments of straight lines, 

b) Segments of relatively large radius circles, 

c) Inflection points, 

30 d) Points of maximum numerical value of the curvature (corners), 

e) Points separating portions of very low and very high numerical value of the curvature, 

f) Small area entities enclosed by a contour. 

A set of two or more primitives with specified characteristics is called a feature, illustrated 
35 in figure 2b. 

As previously mentioned there are some requirements to useful features: 
1) A feature should preferably have a reference point. 
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2) A feature should preferably have a unique direction in the image. 

3) A feature should preferably have one or more rotation invariant descriptors. The 
properties described by such descriptors will be called intrinsic properties. 

5 The requirement 3) is not strict. In the absence of intrinsic properties, the match search 
will be different. In this case the comparison between intrinsic features of training and 
recognition images is cancelled and the recognition is based solely on the cluster search. 

If the Image of an object contains only few features of one kind, additional feature types 
10 should be included in the analysis. 

The aim of this section is to describe and exemplify how primitives and features are 
derived from the contours. 

15 Curvature versus contour length 

A good tool for generating primitives is a function describing the curvature versus arc 
length along the contour. Let the tangent direction at a point on a contour be given by the 
angle 6, and let s be the arc length along the contour measured from an arbitrary 
reference point. Then the curvature is de/ds. The curvature function k(s) = de/ds versus s 
20 is useful for defining primitives. Thus zeros in k(s) and reasonably high values of |dic/ds| 
are inflection points. Positive peaks of k(s) are concave corners, negative peaks of k(s) are 
convex corners (or opposite depending of the definition of background and foreground). 
Straight sections of the contour has k(s) « 0 in a range of s. Circular sections with radius R 
has k(s) = +/-1/R in a range of s. 

25 

The due to the pixel discretization, the functions e(s) and k(s) are derived by replacing 
differentials by differences. For this to be meaningful it is preferred to work with high 
accuracy and efficient noise reduction. Sub-pixel definition of contours is essential (see 
Figure 4), and image blurring is often necessary in order to reduce the camera noise. It is 
30 also helpful to smooth the contour function k(s) before deriving the primitives. 

Figure 5 shows the behaviour of the curvature function k(s) in case of the outer contour of 
the image in figure 3. In figure 6 are shown the symbols for straight sections and corners 
detectable using the curve in figure. 5. There is a similar correspondence between zeros in 
35 k(s) and inflection points (not shown). 

The algorithms for generating primitives need certain threshold values for the curvature. 
For example a straight line is characterized by |k(s)|< K a over a range of s, where k 3 is the 
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curvature threshold, and the integral kds over the range should also be sufficiently small 
(below an angle threshold) since kds represents the tangent angle variation. Another 
threshold K b is relevant for deciding if a positive or negative peak is a corner or just noise. 
The corner criterion is then: [K(s)>K b/ and k(s) is the local maximum] or [k(s)<-k 6/ and 
5 k(s) is the local minimum]. 

Figure 5 illustrates Curvature k(s) (in radians/pixel) as a function of arc length s (in pixels) 
along one of the contours of Figure 3, Window 4. The symbols are used for showing the 
correspondence, see figure 6. Figure 6 illustrates the outer contour found for the image 
10 shown in Figure 3. The symbols used for characteristic features correspond to those given 
in Figure 5. 

Steps in the 3D Interpretation 

In subsection (The training process) below, the training process is described in which a 
large number of training images is created. In subsection (The recognition process) below 

15 the steps in the recognition process are outlined, i.e. how features derived in the training 
process are compared with those of a recognition object. In Section (The recognition 
process in case of pairs of line segments as features) the steps of the recognition process 
in special case of line pairs as features are described. The special status of the distance 
parameter p in the object-camera pose and the use of two or more cameras are discussed 

20 in Section (The special status of the parameter p, the use of two cameras). 

The training process 

In the training process a large number of images of the object with known object-camera 
poses are generated. This can be done by construction in a CAD system or using a camera. 

25 

The training geometry is illustrated in figure 16, wherein the frame of the object is given 
by the axes x,y,z. The optical axis of the camera is going trough the origin of the object 
frame. The horizontal u-axis of the camera is preferably parallel to the x-y-plane of the 
object frame, see figure 16. The training parameters are p,cp,6. 

30 

It is necessary to produce many training images corresponding to different camera poses 
(positions and orientations) relative to the object. Because of the use of: 

1) homographic transformations during recognition and 
35 2) rotation invariant intrinsic descriptors, 
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the training Involves only 3 degrees of freedom. These degrees of freedom are chosen to 
be the spherical coordinates (p,<p,6) of the camera pinhole in the frame of the object. The 
angular pose of the camera is characterized by an optical axis going through the origin of 
the object frame, and a horizontal image axis parallel to a specified plane in the object 
5 frame (see section, 'Recognition using tilt-pan homographic transformation 1 ). The camera 
poses used in the training are suitably distributed in the p,cp,8-space. Usually the chosen 
poses form a regular grid in this space. The discretization step of cp and 6 is preferably of 
the order 2-5 degrees. The range of p and the number of different p-values depend on the 
situation. In this presentation we do not go in detail with the distribution in p,(p,0-space of 
10 the training poses. A single index i is used for the training pose assuming a well-known 
relation between this index and the corresponding pose of the training camera relative to 
the object. 

The flow chart in figure 11, describes the processing for training. 

15 

The recognition process 

Consider a certain feature type. Each training image contains a number of features. Let 
these features be n M where i is the index of the training image and j is the index of the 
feature in the image. Let now «, be feature j in the recognition image. Each feature has the 
20 following properties: 1) the reference point Q, 2) the reference direction defined by the 
angle y, and 3) intrinsic properties consisting of one or more numerical quantities. We 
denote the intrinsic properties by the vector A. Note that y and the components of the 
vector A must be derived after a tilt-pan homographic transformation (see section, 
'Recognition using tilt-pan homographic transformation 1 ) moving the point Q to the mid- 
25 image point. 

The match search Involves a 1) comparison of A(7ij) with A(n,j) and 2) a cluster search in 
the parameter space describing the pose of the potential recognition objects. If A^j) is 
sufficiently similar to A(n u ) then there is match in relation to intrinsic parameters. For all 

30 intrinsic matches, the recognition pose derived from i, Q(n,,j), y(nij) ,Q(nj), y(n$ is 

calculated and used in the cluster search. Here, three degrees of freedom of recognition 
pose is given by the index i defining the training pose, while the other three degrees of 
freedom are tilt, pan, and roll in the relevant homographic transformation. The 
mathematical details of this step are given in Sect. (Recognition using tilt-pan 

35 homographic transformations). 

Each accepted cluster could be considered to represent a physical object. However, 
additional checks for 3D overlap between the guessed poses should be performed after the 
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cluster search. The (p,<p,6) configuration space for training is necessarily discretized, and 
so a simple recognition procedure gives an error of the order one half discretization step. 
This error may be reduced by interpolation between results from neighbour training 
images. 

5 

Figure 17 and 18 illustrates the structure of the database of descriptors derived from 
training images and the structure of descriptors derived from recognition image. Pi,<pj,e k are 
descretized values of training parameters. Each record (line) in the tables is derived from a 
feature. In the present example each feature has three intrinsic descriptors. The number of 

10 extrinsic descriptors is preferably 3: The two coordinates of the reference point of the 
feature and the reference direction after homographic transformation. Any extrinsic 
descriptor of a record of the database and any extrinsic descriptor of a recognition record 
define together a tilt-pan-roll transformation of brining the recognition feature to coincide 
with the training feature. An increment at the corresponding point in tilt-pan-roll 

15 parameter space can then be performed. In case that the intrinsic recognition descriptors 
are sufficiently different from the intrinsic descriptors of the database the corresponding 
pair of features are not considered. This omission reduces the noise from false 
correspondences. 

20 Many important details in this recognition process are presented in the example given in 
the next section. 

The recognition flow chart is shown in figure 12. 

25 The recognition process in case of pairs of line segments as features 

In the following example the features used in the recognition are pairs of line sections, and 
the object is a 'triple brick' consisting of 3 box formed elements, see figure 7. Figure 7 
illustrates the 3D object treated in the example of the pose determination. 

30 Straight sections are preferably derived from the contours of the training photos. All pairs 
of line sections in each training image are then considered as a feature. The intersection 
point of the line pair is the point Q, the angle y Is the angle between the horizontal image 
axis and the direction of the bisectrix of the line pair. Intrinsic descriptor are: 1) the angle 
V between the lines, and 2) the distances between intersection and the end points of the 

35 line section. Both types of intrinsic features are derived after the homographic 
transformation. 
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The distances between the Intersection and end points should not be used directly In the 
match search, since partial occlusion of a line segment produce wrong distances. 

The following discussion is focusing on searches based on the angular descriptor only. 

5 

Figure 8 illustrates training photos of Lego model of the triple brick structure. 

Figure 9 illustrates the curvature in units of radians/pixel as a function of the arc length (in 
pixels) along the contour in figure 8A. 

10 

As seen in figure 9 it is easy to localize straight sections. There are 12 straight sections in 
this example. When deriving relevant pairs of lines some pairs are omitted, namely those 
with angles near 0 or 180 degrees between the lines. In this way a total of about 90 pairs 
in figure 8A can be considered. 

15 

Figure 10 shows a combined view of training and recognition image. P is the focal center. 
The training line pair ab is homographically transformed so that Q ab moves to P. This gives 
rise to the line pair a'b\ The quantity y ab is the angle between the horizontal direction and 
the bisectrix m ab . The intrinsic descriptor of the ab line pair is the angle V ab between a' and 
20 b\ Similar definitions hold for the quantities derived from the line pair cd of the recognition 
image. Intrinsic descriptors other than the angle V are not shown. 

Figure 10, shows a combined view of a training image and a recognition image. In the 
example the positions Q ab and Q cd , and the angles y cd and y ab between the bisectrices m cd 

25 and m ab define a tilt-pan-roll transformation between the training pose and the recognition 
pose. Thus a line pair of the recognition image related to a line pair of a particular training 
image defines a point in tilt-pan-roll parameter space. The line pair ef shown at the top 
right part of figure 10 has nearly the same angular descriptor as the ab line pair. And so, 
in the algorithm this line pair comparison produces a 'false' point in parameter space. 

30 However, it is characteristic for non-corresponding primitives that they produce very 
scattered points in parameter space, while corresponding line pairs gives a cluster in 
parameter space. If another training image is attempted to match with the recognition 
image in figure 10, one does not get any clustering of matching line pairs. 

35 It is clearly seen that object occlusion and insufficient image segmentation does not harm 
the recognition process unless the background of false points in the parameter space 
becomes comparable to the signal formed by clusters of true matches. 
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If the number of false matches disturb the cluster search, it is possible to setup 
inequalities involving the above-mentioned additional intrinsic descriptors (distances 
between end points of line sections and the intersection point). Such inequalities allows the 
recognition lines to be only partially visible, but forbid the recognition lines to have any 
5 section present outside the range of the training lines. 

Given a recognition images with features % s the algorithm runs as follows: 

For each training image index i do 

{ 

10 Reset roll, tilt, pan parameter space; 

For alt valid indices j and j' compare A(n u ) and A(n r ) do 
{ 

In case of intrinsic match 

15 Derive roll r, tilt t, pan p, from Q(n {<j ), y(n y ) ,Q(n y ), Y(ny); 

Update parameter space accordingly; 

> 

> 

Test for clustering and store coordinates of clusters with sufficiently high 
20 density/population along with the index i of the training image. 
} 

Above the "intrinsic match' is based on similarity of the angles V and fulfilment of 
inequalities concerning the distances. The term 'update parameter space' means to 
25 increment the vote at the relevant point in parameter space. 

In case of weak perspective, a single object may produce clusters with several training 
images, namely those having equal cp- and 0-values, but different p-values. Only one of 
these p-values corresponds to a real object, and so a special algorithm using the linear 
30 intrinsic descriptors should be used. 

The back-transformation using the guessed training index i, and the cluster point in the 
tilt-pan-roll space is described in Sect. (Recognition using tilt-pan homographic 
transformations). The recognition process is now completed, since this back- 
35 transformation defines a pose of the recognition object relative to the camera (or visa 
versa). 
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The special status of the parameter p, the use of two cameras 

The object-camera pose fs described by 6 parameters, namely a) p,<p,8 of the 
corresponding training image and b) (roll,tilt,pan) of the transformation between 
recognition pose and training pose. In case of weak perspective, the dimensionless 
5 descriptors of the primitives (such as angles) are almost independent of p, and the linear 
descriptors are approximately proportional to 1/p. Therefore, with weak and moderate 
perspective It is possible to limit the training to rather few different p-values, and use a 
suitable interpolation. Preferably the recognition can be split into two parts. The first part 
concentrates on finding the 5 angular parameters, the second part derives the p-value 
10 using interpolation. 

The accuracy of the final value of p depends on the quality of the linear descriptors of the 
features and the number of different p-values involved in the training. Even using 
interpolation, the relative uncertainty of p is preferably not smaller than the relative 
15 uncertainty of those intrinsic descriptors involving distances between image points. 

Increased accuracy can be obtained by using two (or more) cameras since 3D triangulation 
can be carried out as follows. 

Consider two cameras Camera 1 and Camera 2 (Fig, 13). Let a classification and pose 
20 estimation result of Camera 1 be characterized by (id, p lf cpi,0i, t lt p u r x ) where \c t is the 
object type index and the remaining parameters define the pose. If p ± is completely 
uncertain, the remaining pose parameters define a line L x for the object reference point 
(see Fig 13). Let in a similar way a guess (ic 2 , p 2 , q> 2 , &2/ 1 2 , p 2 , r 2 ) of Camera 2 define a 
line L 2 . Since the global camera poses are known, the lines L 2 and L 2 can be represented in 
25 the global frame. In case that iCi and lc 2 are equal, and the L x and L 2 essentially cross each 
other, then the pair of guesses is assumed to represent a real object. Furthermore, the 
two - previously uncertain - parameters pj and p 2 can be determined by the pseudo 
intersection points (Fig 13) with high accuracy. In this way one can not only enhance the 
accuracy of the pose estimation, but also avoid misinterpreting false single-camera results. 

30 

In the following recognition using tilt-pan homographic transformations will be described in 
greater details. 



WO 2005/017820 



24 



PCT/DK2004/000540 



Recognition using tilt-pan nomographic transformations 



We shall use the following notation: 



10 0 

l? x (/)=<jo cos/ sin? 
0 -sin/ cos/ 



cos p 0 - sin p 
R y (p)=' 0 1 0 
sin p 0 cos p 



cosr 


sin;- 


0' 




f 


0 


0^ 


« -sinr 


cosr 


0 


K = < 


0 


/ 


0 ► 


0 

v. 


0 


1 

J 




0 

v. 


0 


1 

J 



(1) 



5 The R matrices describe rotations through the tilt angles t (tilt), p (pan) and r 
(roll). The matrix K describes the perspective transformation with the focal length 
f. We shall use the composite rotation 5R(/ 5j p,r) defined as 



* (t,P,r)=RAr)R y (p)R x {t) 



(2) 



10 



where the succession of the matrix multiplications is running from right to left. 
^(^•^^(-^ (-^.(-r) is the inverse transformation. 



The elements of $l(t,p,r) are 



15 



cos r cos p sin t sin p cos r - cos t sin r - cos t sin p cos r + sin t sin r 
- sin r cos/? -sin/sin psinr + cos/cosr cos/sin psinr + sinfcosr 
sin /? - cos p sin f cos p cos / 



We assume that a combined rotation-translation of a camera relative to an object 
frame is defined by %l{$ 9 p 9 r\D, where D is the translation vector expressed in the 
20 object frame. Then, if a point is represented by coordinates (x 0 y 0 z 0 ) in the object 
frame and the same point is represented by the coordinates (x c y c z c ) in the camera 
frame, one has 
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fx > 
x c 


= 9*(f ,/?,/) 








\ Z C J 







The coordinate axes of the camera frame expressed in the object frame are row 
vectors of 5R(/, />,;•). Therefore the camera z axis has the global direction 
(sin p,- cos p sin t, cos p cos t) . 

5 

We now define the training pose of the camera. We chose a pinhole position 
given by D = {- p sin q>, p cos <p sin 6,-p cos <p cos 9) . We chose the angles (t, p, r) of 
the training camera to be equal to {p,(p,yQ. This implies that the optical 

axis goes through the object origin. The choice r = im P |ies that the x c 
10 axis of the camera is parallel to the yz-plane of the object frame. 

Let the object-camera transformation be T^T ■ This transformation shall be split 
into two contributions: 



15 



20 



y Camera _ ^Camera 'T'TraimngCamera 
Object ^TrainingCamera 1 Object 

Tar* 5 — = n(e,<p,-%) \D{p,<p,e) 

D(p, <p, 0) = (- p sin <p, pcos<psm0,-p cos <p cos B) 



(3) 



(4) 



where O is the zero vector. The pose of the training camera is characterized by 
1) an optical axis oriented towards the origin of the object frame, and 2) a camera 
x c -axis parallel with the object yz plane. The angles r, t, p are roll, tilt, and pan 
angles characterizing the transformation T^^^from the training camera 
25 orientation to the recognition orientation. Note that the total transformation TST™ 

Object 

can be derived from p and the five angles r, t, p, B, and q>. 
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Let Q.{p,<p,e,r,p,i) be an image recorded by a camera with the pose given by 
(p,<p,0,r,p,t). The relation between homographies with common pinhole position, 
i.e. between, say, Cl{p,<p,0,r,p,t) and n(p,<p,0,O,O,O) is a 2D transformation here 
expressed using homogenous coordinates 



cu 
cv' 
c 



= K3i{t,p,r)K-' 



(5) 



where (u, v) and (u' # v'; are image coordinates. The origin of the image coordinates 
is assumed to be at focal point. 

10 We shall focus on the transformation KSi^p^K-' (zero roll angle). The reason is 
that K and R. commute so the roll can be considered as an image rotation. The 
homographic (2D) operation H(t, p) defined by H(t,p) = Km(t,p,6)K- 1 can be 



written 



15 



(„. v .) = A u cosp + v sin psint — f cos / sin p), f(v cost + f sin t ) 
u sin p - v sin t cos p + f cos t cos p 



(6) 



Note that the operator H(t, p) can operate on image points as well as on a full 
image. A transformation moving a specific point Q = (u Q ,v Q ) to the origin (the focal 
point) is given by pan and tilt angles 



t 0 (Q) = - arctan ^r, Po (Q) = arctan . U ° 
f 



(7) 



20 The inverse transformation is given by: 

( M> v) = /(»' cos P ± f sin p), f(u' sin p sin t - v' cos t ~ f sin t cos p) 

-u'smpcost + v'smt + fcostcosp ' 

Let us consider a feature which in the relevant images has 1) a reference 2D-point 
Q and 2) a direction angle y with the u-axis. This feature has the parameters 
(Q lr ,y lr ) in the training image0, r and {Q rec ,y nc ) in the recogition image<3„ c . Thus 
25 Q lr and Q nc are corresponding points and the directions with angles y )r and y nc are 
also corresponding. It is interesting to compare the two images 



WO 2005/017820 



27 



PCT/DK2004/000540 



These two images are rotated versions of each other, and the rotation angle of Ci' lr relative 
5 t0 ^'rec is Y n c ~7, r - They both correspond to poses with the optical axis going 
through a 3D point Q which is imaged in Q, r , and Q rec . It can be proved that the tilt, 
pan, and roll angles, t, p, and r of the recognition pose relative to the training pose is 
given by: 

io *(r,/v-)= ^(fo^lPoiQ^rJ^ifoiQ^Poi&rlr.) (io) 

The treatment of the training images runs as follows: The images are analysed in 
terms of features (single fixpoints, fixpoint groups, composite features). The values 
Qrec, ynec and various other intrinsic numerical descriptors of each feature are derived 
15 and stored in a database. For speed optimizing purpose, the elements of the matrix 

is also calculated and stored. 

During recognition, the intrinsic numerical descriptors of all features in the 
recognition image is compared with the descriptors of all similar features in the 
20 database. In case of sufficient match the parameters t, p, r of the rotation T<i amera 
are derived from (9). The resulting values of t, p r and the parameters p,<p,9 of 
the matching training image form together a guess of a 6-parameter pose of a 3D 
object, since they define the transformation . Clusters in the 6-dimensional 

pose space are assumed to be caused by a real object located in the corresponding 
25 pose. 

In case of weak or moderate perspective, the number of different values of the 
parameter p may be chosen to be very small, since length features has a p - 
dependence approximately equal to p' 1 . In this case it is recommendable to let the 
30 scale-invariant descriptors determine the relevant point in the 5 dimensional 
r,p,t,q>,B -space. A subsequent analysis of the scale variant descriptors then 
determines the relevant p -value. 
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The total object-camera transformation can be expressed as: 



rn Object 

■* Camera 



n(t, P ,r) \o n[e,<p,y$ \D{ P ,<p,e) 



